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We derive the relation between the amplitudes of short-baseline appearance and disappearance 
oscillations in 3+iV s neutrino mixing schemes which is the origin of the appearance-disappearance 
tension that is found from the analysis of the existing data in any 3+N s neutrino mixing scheme. 
We illustrate the power of the relation to reveal the appearance-disappearance tension in the cases 
of 3+1 and 3+2 mixing using the results of global fits of short-baseline neutrino oscillation data. 
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The standard three-neutrino mixing scheme, which ex¬ 
plains the oscillations that have been observed in so¬ 
lar, reactor and accelerator experiments (see m Ei), 
may be an incomplete description of neutrino mix¬ 
ing, as suggested by the indications in favor of short- 
baseline oscillations generated by a squared-mass dif¬ 
ference Am| BL ^ 1 eV 2 which is much larger than the 
standard “solar” (Am| OL ~ 7.5 x 10 -5 eV 2 ) and “atmo¬ 
spheric” (A m\ TM & 2.4 x 10 -3 eV 2 ) squared-mass dif¬ 
ferences characteristic of three-neutrino mixing. 

The indications in favor of short-baseline oscillations 
are: the reactor antineutrino anomaly [3], which is a 
deficit of the rate of v e observed in several short-baseline 
reactor neutrino experiments in comparison with that ex¬ 
pected from the latest calculation of the reactor neutrino 
fluxes HE]; the Gallium neutrino anomaly mm , con¬ 
sisting in a short-baseline disappearance of v e measured 
in the Gallium radioactive source experiments GALLEX 
[TT1 and SAGE [12]; the signal of short-baseline -+ v e 
oscillations observed in the LSND experiment [13, 14]. 

These results can be explained by extending the stan¬ 
dard framework of three-neutrino mixing with the in¬ 
troduction of new massive neutrinos that in the flavor 
basis correspond to sterile neutrinos [15], which do not 
have standard weak interactions. The simplest scheme is 
that called “3+1”, in which there are the three stan¬ 
dard active neutrinos z+, v T , and one new sterile 

neutrino u s . In the mass basis, there are the three 
standard massive neutrinos z+, z/ 2 , ^3, whose two inde¬ 
pendent squared-mass differences correspond to Am| OL 
and A rn\ TM , and an additional massive neutrino z/4 
which generates the short-baseline squared-mass differ¬ 
ence Acx A m\ 2 — Ara 2 3 cx Am| BL (with the usual 
definition Am^ = ra 2 — ra|, where rrik is the mass of 

^)- 

However, it is possible that there are more than one 
new massive neutrinos which generate more than one 
Am| BL A) leV 2 . In the literature one can find stud¬ 
ies of 3+2 [EH20], 3+3 [H], 3+1+1 [2TH21], and 1+3+1 
[25] [26] schemes. It turns out that in all these schemes 
there is a tension between the results of short-baseline ap¬ 
pearance and disappearance experiments 


S3]. The reason of this tension is well understood in the 
framework of the 3+1 mixing scheme from the analytic 
relation which connects the amplitudes of appearance 
and disappearance short-baseline oscillations [27][28]. In 
this paper we extend this explanation to the general case 
of 3+7V s mixing, in which there are N s mostly sterile 
massive neutrinos at the eV scale (some of which can 
also be at a larger scale). 

For the study of neutrino oscillations in vacuum it is 
convenient to use the following general expression of the 

probability of ( z+* oscillations [451 16] : 

(-) = fiot/3 d ^ |£7/3 /c| ) sin A ]$p 

/c+p 

+ 8 ^ ^ | U a jUpjUakUpk | sin A/^p sin A^p 

j>k 

(+) 

x cos(A jfe - rj a p jk ), (1) 

where 

Am? L 

A/cp = , Ooifijk = ar S tyotjUpjUotkUpj^ , (2) 

and p is an arbitrary fixed index, which can be chosen 
in the most convenient way depending on the case under 
consideration. In the case of three-neutrino mixing, there 
is only one interference term in Eq. 0, because for any 
choice of p there is only one possibility for j and k such 
that j > k. 

In the following we take into account that the non¬ 
standard massive neutrinos must be mostly sterile, i.e. 

\U a k\ 2 < 1 (a = e, n, t; k = 4,...,N), (3) 

in order not to spoil the successful three-neutrino mixing 
explanation of solar, atmospheric and long-baseline neu¬ 
trino oscillation measurements. The non-standard mas¬ 
sive neutrinos are conveniently labeled in order of increas¬ 
ing mass: < ... < mjy. 

We are interested in the effective oscillation proba¬ 
bilities in short-baseline experiments, for which A21 < 
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FIG. 1. Three-dimensional illustration of the relation in 
Eq. § (without approximation) and the approximated rela¬ 
tion in Eq. ( |10| (with approximation). 

A31 <C 1. In the general case of 3 +N S mixing with 
Am|i > Ara| BL an d ^ki > 1 for & > 4, choosing p = 1 
in Eq. Q we obtain 

N 

Ff F - S aP - 4^ |t/«+ (<S Q/J - |C/+ 2 ) sin 2 A*,! 

Va-*vp k=4 

N N 

+ 8 E E I U a jUpjUakUpk | sin Afei sin A71 
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x eos(Aj/ c - rj a/ 3 jk ). (4) 
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FIG. 2. Allowed regions in the sin 2 2$^-sin 2 2#j^ plane ob¬ 
tained in the total 3+1 fit of short-baseline data correspond¬ 
ing to the 3+1-LOW fit of Ref. |33]. The green shadowed 
regions are the regions allowed at lcr, 2cr and 3cr by the anal¬ 
ysis of short-baseline disappearance (DIS) data, with the best 
fit value indicated by a dark-green cross. The strips enclosed 
by the blue diagonal lines are allowed at lcr, 2cr and 3cr by 
the analysis of short-baseline appearance (APP) data, with 
the central best fit dark-blue line. The solid lines correspond 
to the exact relation in Eq. ^ , whereas the das hed lines cor¬ 
respond to the approximated relation in Eq. (10). The regions 
inside the red-orange closed curves are allowed at lcr, 2cr and 
3cr by the global (GLO) analysis of short-baseline data, with 
the best fit value indicated by a dark-red cross. 


transition amplitudes 


Considering the survival probabilities of active neutri¬ 
nos, let us define the effective amplitudes 

sin 2 2+2 = 4| U ak \ 2 (1 - | U ak \ 2 ) =* 4|C/ a fe| 2 , (5) 

for a = e, /i, r and k > 4. The approximation is due to 
the constraint ©, which allows to neglect the quadrat- 
ically suppressed contribution proportional to \U a k\ A - 
Dropping the quadratically suppressed terms also in the 
survival probabilities, we obtain 


N 


- 1 - E sin2 2 +« sin2 Afci > 


(6) 


fc=4 


sin 2 2+2 = MUak\ 2 \U/3k\ 2 , (7) 

for t - /3 and k > 4, which allow us to write the transi¬ 

tion probabilities as 

N 

+)+) - E sin2 sin2 Afci 

Votin'J3 k =4 

N N 

+ 2^7 E sin 2 ^S sin 2 ^a/3 sin Afel sin A?' 1 

k= 4 j=k-1-1 

(+) 

x cos(A jfe - T] a p jk ). (8) 


for a = e,n,T. Hence, each effective mixing angle +2 
parameterizes the disappearance of u a due to its mixing 
with ( z+ for k > 4. 

Let us now consider the probabilities of short-baseline 

is a vp transitions between two different active neutri¬ 
nos or an active and a sterile neutrino. We define the 


From the first line one can see that each effective mixing 
angle parameterizes the amount of ( z+, -+ ( z/^ transi¬ 
tions due to the mixing of ( z+, and ( z+ with for k > 4. 
The second line in Eq. (J8| is the interference between the 

v k and Vj contributions for k,j > 4, which depends on 
the same effective mixing angles. 
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FIG. 3. Marginal Ay 2 as a function of sin 2 obtained in 
the total (upper panel) and pragmatic (lower panel) 3+1 anal¬ 
yses of short-baseline data presented in Ref. [44]. The hori¬ 
zontal lines gives the values of A% 2 corresponding to 68.27% 
C.L. (lcr), 90% C.L., 95.45% C.L. (2a), 99% C.L. and 99.73% 
C.L. (3a) for one degree of freedom. The curves obtained 
from disappearance (DIS) data are equal in the two panels. 
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FIG. 4. Allowed regions in the sin 2 sin 2 2f+ plane 

obtained in the pragmatic 3+1 fit of short-baseline data cor¬ 
responding to the 3+1-HIG fit of Ref. [44]. See the caption 
of Fig. [2] The green shadowed regions obtained from disap¬ 
pearance (DIS) data are the same of Fig. [ 2 ] 


FIG. 5. Marginal Ay 2 as a function of sin 2 2obtained in 
the total (upper panel) and pragmatic (lower panel) 3+2 fits 
of short-baseline data presented in Ref. m See the caption 
of Fig. [ 3 ] 

Considering now the transitions between two different 
active neutrinos, from Eqs. © and 0 one can see that 
for each value of k > 4 the transition amplitude sin 2$^ 
and the disappearance amplitudes sin 2$22 and sin 2$^ 
depend only on the elements in k th column of the mixing 
matrix and are related by 

sin 2 20<§ = ^1 - yj 1 - sin 2 2+^ 

x { 1 ~ y / 1 -^ 2 + 3 ^ > ( 9 ) 

for a = e, /i, r. Taking into account the constraint ©, 
we have 

sin 2 20$ ~ 1 sin 2 20« sin 2 20$, (10) 

as illustrated in Fig. [l] This relation was derived in the 
case of 3+1 mixing and k = 4 in Refs. nang. Now we 
see that the same relation is valid in the general case of 
3 +N S mixing between the amplitudes of appearance and 
disappearance oscillations due to each squared-mass dif¬ 
ference Am^ generated by a non-standard mostly sterile 
massive neutrino at the eV scale (or at a larger scale). 
Note that the relation is very powerful in the global fit of 
short-baseline oscillation data, because it constrain the 
appearance and disappearance amplitudes for any fixed 
value of the corresponding squared-mass difference. 

The relation © is very important, because it con¬ 
straints the oscillation signals that can be observed 
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FIG. 6. Marginal Ay 2 as a function of sin 2 obtained in 
the total (upper panel) and pragmatic (lower panel) 3+2 fits 
of short-baseline data presented in Ref. m- See the caption 
of Fig. [3] 


in short-baseline appearance and disappearance experi¬ 
ments in any 3+N s mixing scheme with sterile neutrinos. 
Its experimental test is crucial for the acceptance or rejec¬ 
tion of these schemes. In particular, since both sin 2 2$^) 
and sin 2 are small, the amplitudes of short-baseline 

^ ( ^e transitions implied by the results of disap¬ 
pearance experiments are quadratically suppressed and 
in tension with the short-baseline -+ u e transitions 
observed in the LSND experiment mmmmm- 
In the following we illustrate the importance of the con¬ 
straint © (and its approximation ©) in the cases of 
3+1 and 3+2 mixing using the results of the fits of short- 
baseline neutrino oscillation data presented in Ref. [43] . 
We denote with the labels “APP” and “DIS” the fits of 
appearance and disappearance data alone, respectively. 
The label “GLO” denotes the global combined fit of ap¬ 
pearance and disappearance data. Furthermore, we dis¬ 
tinguish between “total” and “pragmatic” fits. A total 
fit corresponds to a “LOW” fit in Ref. [44], in which 
all the data of short-baseline neutrino oscillation experi¬ 
ments considered in Ref. [44| are taken into account. A 
pragmatic fit corresponds to a “F1IG” fit in Ref. [441, in 
which all the data of short-baseline neutrino oscillation 
experiments considered in Ref. [44] are taken into ac¬ 
count except the anomalous low-energy bins of the Mini- 
B 00 NE experiment SZl EHl. This “pragmatic approach” 
was advocated in Ref. [44] because the anomalous low- 
energy bins of the MiniBooNE experiment f47, [48] are 
incompatible with neutrino oscillations and may be due 
to background (this problem is going to be investigated 


FIG. 7. Allowed regions in the sin 2 2$i^-sin 2 2$i^ plane ob¬ 
tained in the total 3+2 fit of short-baseline data correspond¬ 
ing to the 3+2-LOW fit of Ref. [33]. The green shadowed 
regions are the regions allowed at lcr, 2cr and 3 a by the anal¬ 
ysis of short-baseline disappearance (DIS) data, with the best 
fit value at sin 2 2$i 4 J = sin 2 2$^ = 0. The regions inside the 
blue closed curves are allowed at lcr, 2cr and 3cr by the analy¬ 
sis of short-baseline appearance (APP) data, with the best fit 
value indicated by a dark-blue cross. The regions inside the 
red-orange closed curves are allowed at lcr, 2cr and 3cr by the 
global (GLO) analysis of short-baseline data, with the best 
fit value indicated by a dark-red cross. 


in the MicroBooNE experiment 03 ED])- Note that only 
the APP and GLO fits are different in the total and a 
pragmatic fits, whereas the DIS fit of disappearance data 
is the same. 

Figure [ 2 ] shows the regions in the sin 2 2'#e?-sin 2 2$^ 
plane allowed at lcr, 2cr and 3cr by the total 3+1 fit of 
short-baseline data presented in Ref. [33] (corresponding 
to the “3+1-LOW” fit in Ref. [33] ). One can see that 
there is a clear tension between the regions allowed by 
disappearance and appearance data. The relation § 
allows us to illustrate the tension also by showing the 
marginal Ay 2 as a function of sin 2 2$^ in the upper 
panel of Fig. [3] One can see that there is no overlap 
of the intervals of sin 2 2allowed by the appearance 
and disappearance data at less than about 2cr and the 
result of the global fit is a compromise with the best-fit 
dominated by the disappearance data. 

The fit of the anomalous MiniBooNE low-energy bins 
with 3+1 neutrino oscillations requires a small value 
of A mh and a large value of sin 2 2$i 4 J , which are 
in strong tension with the data of disappearance ex¬ 
periments 03 S3 m. Hence, we can alleviate the 
appearance-disappearance tension by adopting the prag- 
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FIG. 8. Allowed regions in the sin 2 2$i)/-sin 2 2$^/ plane 
obtained in the pragmatic 3+2 fit of short-baseline data cor¬ 
responding to the 3+2-HIG fit of Ref. [33]. See the caption 
of Fig. [7| The green shadowed regions obtained from disap¬ 
pearance (DIS) data are the same of Fig. [t| 


matic approach described above. The allowed regions in 
the sin 2 2'#e?-sin 2 plane in Fig. [I] obtained with the 
pragmatic 3+1 fit of short-baseline data show that the 
appearance-disappearance tension is acceptable. A com¬ 
parison with Fig. [2] shows that the strips allowed by the 
analysis of short-baseline appearance (APP) data have 
moved towards smaller values of sin 2 2^) and sin 2 2 $j^, 
which are in agreement with the analysis of disappear¬ 
ance data. The same behavior can be observed by con¬ 
fronting the two panels in Fig. [3} where one can see that 
in the pragmatic approach the intervals of sin 2 2 $^ al¬ 
lowed by the appearance and disappearance data overlap 
at about Icr. 

Let us now consider the 3+2 mixing scheme. Figures [5] 
and ^ show the marginal Ay 2 as a function of sin 2 2$+* 
and sin 2 2^J , respectively. They are different because 
the ordering Am^ < Am § 1 breaks the degeneracy be¬ 
tween the mixing of z /4 and z+. Thanks to the relation <§, 
these figures are sufficient to visualize the appearance- 
disappearance tension. 

The upper panel in Fig. [5] shows that in the total 3+2 
fit there is no overlap of the intervals of sin 2 2$^ allowed 
by the appearance and disappearance data at more than 
2cr, which corresponds to an appearance-disappearance 
tension which is worse than that in the total 3+1 fit il¬ 
lustrated in the upper panel of Fig. [ 3 ] (see the discus¬ 
sion in Ref.|43j)- On the other hand, the upper panel 


in Fig. [6] shows that the data do not constrain much the 
value of sin 2 2$^ and the corresponding marginal A+ 2, s 
do not show an appearance-disappearance tension. How¬ 
ever, it is clear that the manifestation of an appearance- 
disappearance tension in one of the oscillation amplitudes 
is enough. 

The secondary local minimum at sin 2 ~ 1.2 x 

10 -3 of the appearance marginal Ay 2 in the upper panel 
in Fig. [5] is due to the possibility to fit the appearance 
data with sin 2 2$i 4 J zs sin 2 2$^} « 1 — 2 x 10 -3 (with 
A mh « 2eV 2 and Am^ « 4eV 2 ). This can be seen 
from the correlated allowed regions in the sin 2 2‘d^jl- 
sin 2 2$® plane shown in Fig. H in which the 2cr region 
allowed by appearance data has a bulge for sin 2 2'd^l « 
sin 2 2^eJ ~ 1 — 2 x 10 -3 (this type of figure has been 
already presented in Ref. [38]; similar figures for defi¬ 
nite values of Am 2 x and Am^ have been presented in 
Ref. [26]). From this figure one can have another view 
of the appearance-disappearance tension, which is clearly 
dominated by sin 2 2$i 4 J . 

The lower panels of Figs. [5] and [6] and Fig. [8] illus¬ 
trate the appearance-disappearance tension in the prag¬ 
matic 3+2 fit. One can clearly see that, as in the 3+1 
case, it is much milder than that obtained in the total 
fit, indicating an acceptable fit of the data. The sec¬ 
ondary local minimum at sin 2 2$i 4 J « 1.2 x 10 _3 of the 
appearance marginal Ay 2 in the lower panel in Fig. [ 5 ] 
is due to the same reasons as that in the upper panel 
explained above and corresponds to the swelling towards 
sin 2 2$^ ~ sin 2 2$^ ^ 1-2 x 10 -3 of the la region 
allowed by appearance data in Fig. [8] 

In conclusion, in this paper we have derived the rela¬ 
tion © (and its approximation ( [To] )) between the am¬ 
plitudes of short-baseline appearance and disappearance 
oscillations in 3 +N S neutrino mixing schemes. This re¬ 
lation is the origin of the appearance-disappearance ten¬ 
sion that is found from the analysis of the existing data in 
any 3 +N S neutrino mixing scheme. We have illustrated 
the power of the relation to reveal the appearance- 
disappearance tension in the cases of 3+1 and 3+2 mix¬ 
ing using the results of the fits of short-baseline neutrino 
oscillation data presented in Ref. 144- 
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